Charge Density Wave Driven Ferromagnetism in the Periodic Anderson Model 
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We demonstrate the existence of ferromagnetism in the Periodic Anderson Model (PAM) at 
conduction-band filhng near a quarter. We show that this ferromagnetism is not supported by 
Ruderman-Kittel-Kasuya-Yosida (RKKY) interactions but is instead driven by the precursors of 
charge density wave (CDW) formation in the conduction electron band. To study the effect of 
spatial correlations, we compare Dynamical Mean field Approximation (DMFA) and Dynamical 
Cluster Approximation (DCA) results. We find that both RKKY and CDW driven ferromagnetism 
persist as short-range correlations are incorporated into the theory. Both DMFA and DCA show the 
precursors of CDW formation through the strong enhancement of the d-electron CDW susceptibility 
as the temperature decreases, up to the ferromagnetic transition temperature. In addition, the DCA 
captures the signal of a band gap opening due to Peierls instability. 

PACS numbers: 71.27.+a, 75.20.En, 75.30.Mb 
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Rare-earth and actinide-based heavy fermion (HF) sys- 
tems display a wide variety of interesting phenomena, in- 
cluding Kondo insulator, fermi and non-fermi liquid para- 
magnetism, ferromagnetism, antiferromagnetism and su- 
perconductivity : lij . Among these phenomena, ferromag- 
netism is especially interesting since while a number of 
HF ferromagnetic compounds exist ^ , few experimental 
studies on the ferromagnetic phase diagrams are available 
g. Although a number of theoretical studies of HF fer- 
romagnetism are available 0, i, i, 0, i, S [IS El, El, El , 
the range of ferromagnetism in the local moment regime 
has not been fully clarified, and the mechanisms leading 
to ferromagnetic order are still not well understood. 

In this Letter, we investigate charge density wave 
(CDW) driven ferromagnetism occurring near a quan- 
tum critical point (QCP) around quarter filling of the 
conduction band in the Kondo regime of the Peri- 
odic Anderson Model (PAM). The ferromagnetism in 
the region of low conduction band filling (e.g. Ud <C 
0.5) is attributed to the Ruderman-Kittel-Kasuya-Yosida 
(RKKY) exchange interaction between the /-electron lo- 
cal moments mediated by the conduction (c?-) electrons 
0, i, i, 0, i, 1 E3, EH, El, El- Whereas at n<j « 1 
the RKKY and superexchange interactions lead to an- 
tiferromagnetic order Il5| |. However, in the region 
near quarter filling rid ~ 0.5, close to a QCP associated 
with a vanishing ferromagnetic transition temperature, 
the mechanism leading to magnetic order is not well un- 
derstood. We provide theoretical evidence arising from 
the incorporation of spatial correlations into the theory 
that the ordering near quarter filling of the conduction 
band in the PAM is enhanced by CDW fluctuations. Pre- 
vious studies (e.g. and [ll|) likely did not find ferro- 
magnetic ordering at this filling due to the difficulty of 
capturing strong CDW correlations with earlier methods. 

Our interest is enhanced by the discovery of the fer- 
romagnetic superconductor UGe2 where supercon- 
ductivity coexists with ferromagnetism near its ferromag- 



netic QCP. Pressure-dependence studies indicate that 
there may be different mechanisms driving the ferromag- 
netic order in this material [13, [3 • One scenario sug- 
gests that the superconducting-ferromagnetic coexistence 
is related to coupled charge and spin density waves [loj . 
Although our model is not detailed enough to describe 
Uranium-based compounds, our results suggest that a 
CDW mechanism might play an important role in driv- 
ing the ferromagnetic order near the QCP in UGe2. 

As briefly mentioned in a previous study of the infi- 
nite dimensional PAM [3], Kondo screening plays an 
important role in triggering the ferromagnetism when 
the RKKY and superexchange interactions favor anti- 
ferromagnetism. We refer to this as charge density 
wave (CDW) driven ferromagnetism, since the alignment 
of the /-electron local moments is caused by coherent 
Kondo screening induced by the conduction electrons 
forming a CDW. Fig. [1] illustrates this process. Imag- 
ine there are roughly N/2 localized (/) electrons and 
N/A itinerant [d) electrons in the system, where N is 
the number of lattice sites. At some low temperature 
each d electron locally screens one / electron through 
antiferromagnetic alignment of their spins. These local 
screening clouds [H, [2l| become larger at lower temper- 
atures. As a result, the 7V/4 d electrons tend to screen 
all N/2 f electrons coherently. This process requires that 
the d electrons form a spin-polarized CDW, since in order 
to maximize kinetic energy gain and maintain the coher- 
ent Kondo screening a d electron must hop to nearest 
neighbor sites with unoccupied d orbitals. 

Formalism. Since Dynamical Mean Field Theory 
(DMFT) [22*1 in infinite dimensions has a local self en- 
ergy, questions arise as to whether the CDW ferromag- 
netic mechanism persists in a finite dimensional sys- 
tem and how spatial correlations may affect the re- 
sults. In this Letter, we present our study on the three- 
dimensional PAM using the Dynamical Cluster Approx- 
imation (DCA) [23]. We find that the CDW driven fer- 
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FIG. 1: Illustration of a d-electron charge density wave for 
Hf = 1 and Hd — 0.5 which mediates alignment of the /- 
electron local moments through coherent Kondo screening. 
Smaller shaded circles represent the / orbitals, and bigger 
dashed circles the d orbitals. Up and down arrows represent 
the electron spin orientations. 

romagnetism persists as short-range spatial correlations 
are incorporated into the theory. In addition, we present 
one and two-particle evidence of the precursors of CDW 
formation to support this ferromagnetic mechanism. 
The PAM Hamiltonian can be written as 
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H = -t (dld.a + h.c.) + J2i^ddld.a + efflU) 

<ij>(T ia 



(1) 



where d{f)\^J destroys (creates) a d{f) electron at site i 
with spin a, t is the nearest-neighbor d-electron hopping, 
Cd and e f the orbital energy of the d and / electrons re- 
spectively, V the d-J hybridization, and U the on-site 
Coulomb repulsion between / electrons. In this Hamil- 
tonian we set the chemical potential to zero. 

We work with a simple cubic lattice with bare con- 
duction band-width W ^ I2t = i {t ^ 0.25). We 
choose U = 1.5W and V — W/3 to satisfy the Zlatic- 
Horvatic criterion for the strong-coupling regime [2^ : 
u = U/nT > 2.0, where T = V^N{0) and A^(0) is the bare 
conduction band density of states at the Fermi level. The 
Hirsch-Fye Quantum Monte Carlo (QMC) algorithm ^ 
is used to solve the DC A cluster problcm[2f|. We per- 
form calculations for two different cluster sizes: Nc = 1, 
equivalent to DMFA, and Nc — 14. We choose the bi- 
partite cluster 14A, which is considered the first "good" 
cluster with A'c > 1, following the selection criterion of 
Betts [13, [11]. Following the procedure in Ref. for 
every rid value, we choose ed — ef such that n/ « 1 at a 
moderate temperature, say /3 = 1/T — 10, and keep the 
value of e^ — £/ fixed as the temperature is decreased. The 
values of ed and e/ are adjusted (keeping ed — Cf fixed) 
during the self-consistent iterations to satisfy w 1 and 
the chosen value of n^. To study the systematic errors 
associated with the imaginary time increment (At) in 
the QMC we compare the results with tAr — 1/4 and 
tAr = 1/16. 

Results. The upper panel of Fig. [5] shows the FM 
Curie temperature, Tc, versus n^, with Tc vanishing at 
Ud slightly higher than 0.6. Comparison between results 
for Nc = 1 with tAr = 1/4 and tAr = 1/16 suggests 
that tAr = 1/4 is sufficient to capture the correct pro- 
file of the phase diagram. Comparison between results 



FIG. 2: Upper graph: Phase diagram showing Tc as a func- 
tion of rid at n/ = 1. The Tc values are extracted from the /- 
electron ferromagnetic susceptibility data. Lower graph: The 
corresponding AJrkky vs Ud- The ferromagnetic phase dia- 
gram extends into the region where the RKKY coupling is an- 
tiferromagnetic. The ferromagnetic region of RKKY coupling 
extends toward higher Ud values as Nc is increased, suggesting 
that the Fermi surface may be altered due to the introduction 
of k-dependence in the self energy. 



for Nc = 1 and Nc = 14 with the same tAr = 1/4 
indicates that spatial correlations decrease Tc, as gener- 
ally expected [23]. The lower panel of Fig. [2] shows the 
dependence of AJrkky on n^. Here AJjikky is de- 
fined as AJrkky = Jrkky ~ "^rkky^ with Jr^ky — 
jRKKYiq = (tt, TT, tt)) and Jrkky = Jrkky{ci = 
(0,0,0)), so that a negative (positive) AJrkky sig- 
nals a ferromagnetic (antiferromagnetic) RKKY cou- 
pling between neighboring local moments. We define 
jRKKYiq) ~ -2j2^riEk„G'^(k,zc^„)G''(k + q,zc^„), 



with Jfd ~ 8VyU, as in Ref. [15|. For both Nc = 1 
and Nc = 14, we take tAr — 1/4, and show that for 
three different temperatures close to Tc (T=l/60, 1/70, 
and 1/80) there is no significant difference in AJrkky- 
It is clear that for both values of Nc the RKKY cou- 
pling changes from ferromagnetic to antiferromagnetic 
in a region where the phase is still ferromagnetic. In 
this region, where AJrkky > 0, ferromagnetism can 
not be explained by the RKKY mechanism. Note that 
the region of ferromagnetism supported by the RKKY 
mechanism becomes wider as Nc changed from 1 to 14. 
Since, for a given n^, RKKY coupling favors ferromag- 
netic alignment when the magnitude of the Fermi vector 
satisfies kp < 7r/4, this result suggests that the Fermi 
surface shrinks as spatial correlations are introduced. 

To address the mechanism of ferromagnetic order in 
the region of the phase diagram where RKKY interac- 
tions favor antiferromagnetism, we provide evidence that 
in this region ferromagnetism is driven by the precur- 
sors of charge density wave formation. In support of this 
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FIG. 3: The main panel shows the d-electron charge density 
wave (CDW) susceptibihty and the bulk charge susceptibility 
as functions of Ud &t T = 1/60, temperature at which the 
CDW susceptibility is strongly enhanced as indicated in the 
inset. The region of Ud values where the CDW susceptibility 
is greater than the bulk charge susceptibility covers the non- 
RKKY region in the ferromagnetic phase diagram. 



idea, Fig.[3]shows that at low temperatures the d-electron 
CDW susceptibility, xcdw = Xchargeiq. = (7r,7r, tt)), is 
strongly enhanced relative to the 'bulk' (uniform) charge 
susceptibility, XBuik = Xchargei^ = (0,0,0)), for Ud val- 
ues corresponding to the ferromagnetic phase region not 
supported by the RKKY mechanism. Based on geome- 
try one expects the CDW configuration be best achieved 
when Ud = 0.5, which implies that xcdw would be most 
strongly enhanced for Ud = 0.5. However, our results 
show that the maximum enhancement occurs around 
Ud ~ 0.55, which may, again, indicate shrinking of the 
Fermi surface. In the inset of Fig. Owe plot xcdw ver- 
sus T for Hd = 0.4, 0.5, and 0.6, showing its enhancement 
at low temperatures. In these calculations, tAr = 1/4. 
For smaller values of tAr (not shown), xcdw increases, 
but the peak location does not shift. 

To elaborate further on how the CDW driven fer- 
romagnetism works, note that, as shown in the inset 
of Fig. [21 Xcdw reaches a maximum at temperatures 
around Tw 1/40 — 1/60, and decreases slightly when 
the temperature is lowered close to Tc- The non-local 
correlations enhance xcdWj and we also find that they 
enhance the screened local moment for ~ 0.5, while 
they suppress it at other fillings (not shown). This may 
be interpreted as an attempt to balance the kinetic and 
potential energies of the Kondo lattice. The system can 
optimize the exchange energy by localizing a d electron in 
a singlet with the / moment on half the sites. This is bal- 
anced by the kinetic energy in the usual way by allowing 
the d electrons to delocalize and screen other /-electron 
local moments forming Kondo "clouds". However, when 
Ud ~ 0.5, the system may gain additional kinetic energy 
though the gap that accompanies a d-electron charge den- 
sity wave (CDW). Note that this is not a well-defined 
CDW phase, since such an ordering would be manifest in 



the divergence of xcdw, and a full gap would suppress 
Kondo screening. As the temperature is lowered further, 
the local Kondo screening clouds grow, so that the over- 
lap between the Kondo clouds drive the d electrons to 
collaborate to screen the / electrons coherently. In do- 
ing so, each d electron forms a Kondo singlet with an / 
electron at a site only momentarily, then hops to one of 
its neighboring sites to regain the resonance with the / 
electron at that site. This process forces the /-electron 
local moments to align, hence the system forms a ferro- 
magnetic phase. The momentary breaking of the Kondo 
singlet due to hopping of the d electrons to their neigh- 
boring sites decreases their staggered charge correlation. 
Thus, this explains why xcdw is suppressed slightly as 
the temperature approaches Tc- 

So far we have shown that DMFA can already cap- 
ture the d-electron CDW precursors at the two-particle 
level, i.e. in the xcdw behavior. In the following, we 
show that by incorporating the spatial correlations within 
DC A {Nc > 1), the CDW precursors can also be cap- 
tured at the one-particle level, e.g. through the d-electron 
density of states (DOS). In order to capture the effec- 
tive modulation potential inducing the Peierls instability 
[29II, the self energy must have sufficiently strong k de- 
pendence. This is only possible for Nc ^ 1. When this 
modulation potential is properly captured by the self en- 
ergy, the corresponding d-electron DOS forms a gap or 
signals of a gap opening, such as a "dip", around the 
chemical potential. In our model, the chemical poten- 
tial is set fixed to zero, but we adjust to satisfy the 
desired filling. Therefore, if a gap or a "dip" occurs in 
the DOS, it should be manifest in a plot of versus rid- 
Fig. [3] shows how the ed vs rid profile evolves from high 
to low temperatures for Nc — 1 and Nc = 14. It is clear 
that for Nc = 14 (right panel), but not for Nc ^ 1 (left 
panel), the vs curve bends slightly near rid ~ 0.5 
- 0.55 at r « 1/40, which coincides with the strongest 
enhancement in xcdw- We interpret this as the signal of 
a gap opening accompanying the formation of a CDW in 
the d-electron DOS. As the temperature is lowered fur- 
ther, e.g. from T « 1/40 to T « 1/70, the bending in 
the ed vs rid curve does not become more prominent, and 
diminishes as the temperature becomes closer to Tc, e.g. 
T « 1 /90 (not shown) . This indicates that the gap never 
fully opens, and the d-electron DOS increases back as the 
d electrons become more mobile to screen coherently at 
temperatures closer to Tc- 

Our arguments may be justified also by inspection of 
the Kondo screening length [l^, [2l| and its consistency 
with the Xcdw and the d-electron density of states as 
discussed before. The magnitude of the Kondo coupling 
in our system is Jfd « 8V'^/U « 0.6V1^ ~ 7t. According 
to the scaling theory of Sorensen et al. [l^], this cou- 
pling strength roughly corresponds to the Kondo screen- 
ing length, ^L, of the order of a lattice constant or less. 
This is consistent with our picture of localized Kondo 
'clouds' when the CDW forms. This screening length de- 
pends on the d-electron density of states at the chemical 
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FIG. 4: Comparing versus rid for Nc = 14 and Nc = 1. 
For A^c = 14 (right panel) bending in the ea vs rid curve is 
observed, indicating the tendency to the formation of a band 
gap due to Peierls' instability. This feature is not seen for 
Nc = l (left panel). 



potential, p(0), through oc exp[— l/p(0) J/d]- Since 
p{0) iiic;r(;asc;s as the; temperature approaches Tc (as dis- 
cussed in the previous paragraph), the screening length 
also increases accordingly. In turn, the overlap between 



screening clouds enables the d electrons to mediate fer- 
romagnetic coupling between /-electron local moments. 

Conclusion. We have shown that the ferromagnetism 
in the strong coupling regime of the PAM has two mecha- 
nisms. In the region of low conduction band filling ferro- 
magnetism is driven by the RKKY exchange interaction, 
while in the region of higher conduction band fillings up 
to slightly more than a quarter the precursors of a CDW 
formation in the d electrons induce the magnetic order. 
We have demonstrated that the CDW formation can be 
captured within the DMFA and DCA through the en- 
hancement of d-electron CDW susceptibility. In addition, 
the DCA captures the signal of a band gap opening due 
to the Peierls instability that drives the CDW formation. 
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